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O ! ABSTRACT 

(N ■ 

q ' Aims. We have detected a strong limb polarization for Uranus and Neptune. With spectropolarimetric observations we characterize 

, the spectral dependence of this limb polarization and explore the diagnostic potential for investigating the distribution and properties 

/"N . of the scattering particles. 

Methods. We present disk resolved spectropolarimetry of Uranus and Neptune covering the wavelength range from 5300 A to 9300 A 
and compare the spectropolarimetric signal for different limb sections and the center of the planetary disk. As an additional benefit 
we obtained center-to-limb disk profiles for the intensity and polarization for various wavelengths. 
7—i • Results. Our data show for both planets a strong linear limb polarization oriented perpendicular to the limb. The polarization at the 

' center of the planetary disks is essentially zero. Overall the limb polarization increases towards shorter wavelengths, and it is further 

enhanced in the methane absorption bands when compared to adjacent inter band regions. Qualitatively, the polarization signal is very 
t-H , similar for Uranus and Neptune. 

For Uranus the peak limb polarization of the methane absorption bands reaches a maximum of 3 % whereas in the nearby higher 
albedo regions it peaks at about 2 %. The measured polarization for Neptune is significantly lower because the signal is averaged 
t-H ' down by the limited spatial resolution of our Neptune data. 

, The disk profiles for Uranus (center to limb profiles) show for the intensity a strong asymmetry between the bright South pole 

regions and the darker northern latitudes. The asymmetry is particularly strong in the methane bands. We measure a significant limb 
***>; brightening in the methane bands at the Uranian eastern and western limbs. The intensity asymmetry between North and South and 

'"i^ | the limb brightening display a tight anti-correlation with the albedo of the selected wavelength interval. 

From the polarization profiles of Uranus we determine slit integrated polarization values in order to derive intrinsic limb polarization 
parameters which are corrected for the limited spatial resolution of our data. These parameters are compared with previous imaging 
' polarimetry. 

Conclusions. The polarimetric measurements of Uranus and Neptune are compatible with expectations for the limb polarization of 
scattering atmospheres containing substantial layers of particles which scatter radiation similar to Rayleigh scattering. It seems that 
the limb polarization contains interesting diagnostic information, in particular on the scattering properties of the aerosol particles in 
the haze layers at high altitudes. 
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1. Introduction In this second paper we present disk resolved spectropo- 
larimetry of Uranus and Neptune. Our data cover a wavelength 

In a previous paper we have reported the detection of limb po- range from 5300 a to 9300 £ which includes strong methane 

larization of about 1 % for Uranus and Neptune with imaging ab tion bands St absorption bands correspond to strong 

polarimetry in the R, i and z bands (Schmid et al. [20061 here- d jumpfj in the scattering atm0S p her es and corresponding 

after referred as paper i). The linear limb polarization was found polarization effects may be expected for the reflected light, 
to be perpendicular to the limb and roughly equally strong along 
the entire limb. This is unlike Jupiter where the strength of the 
limb polarization is characterized by a strong difference between 
polar and equatorial limb regions (Hall & Riley 1976, Smith & 
Tomasko [19541 Braak et al. 120021 Joos et al. 125031 . 

The limb polarization can be explained as a second order In Sect. [2] we discuss the observations and the data reduc- 

scattering effect which is expected to occur for reflecting atmo- tion. The overall polarization properties of Uranus are described 

spheres seen near phase angle a = 0° (backscattering) where i n Sect. [3] and in the following section we focus on the disk 

Rayleigh-type scattering processes are dominant (see van de profiles (center to limb profiles) for the intensity and the po- 

Hulst|1980l>. larization in some specific wavelength intervals. The spectropo- 

larimetric properties of Neptune are described in Section |3J In 

Send offprint requests to: Franco Joos, e-mail: Sect. the spectropolarimetric data from this paper are com- 

fjoos@astro.phys.ethz.ch pared to the imaging polarimetry presented in paper i and the 

* Based on observations obtained at the ESO 3.6m Telescope at La results are summarized and discussed in the broader context of 

Silla, Chile (ESO program 72.C-0498). scattering polarization from planetary atmospheres. 
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2. Observations and data reduction 

Spectropolarimetric observations of Uranus and Neptune were 
taken on November 29, 2003 with EFOSC2 at the ESO 3.6m 
telescope at La Silla. During the same run also imaging po- 
larimetry of Uranus and Neptune was obtained with the same 
instrument, but with a different setup (see paper i). 

EFOSC2 is a multi-mode imager and grism spectrograph 
which can be equipped with a Wollaston prism and a rotatable 
superachromatic half-wave plate for linear polarimetry and spec- 
tropolarimetry. The Wollaston prism produces together with a 
grism two spectra with opposite linear polarization /y and I ± (the 
ordinary and extraordinary beam) offset by 20" on the CCD. In 
order to avoid confusion a special slit mask is introduced at the 
telescope focus with a series of aligned slits having the lengths of 
19.7" and periods of 42.2". For both planets we have taken data 
with the slit centered on the planet along the N-S direction in the 
celestial coordinate system. For Uranus we have also obtained 
data with slits in the E-W direction which was achieved by ro- 
tating the whole instrument. The spectrographic slits are much 
longer than the diameters of the planets which are about 3.5" 
for Uranus and 2.2" for Neptune. Thus, we obtained "long-slit" 
spectropolarimetry for both planets. The chosen slit width was 
0.5" for Uranus and 1.5" for Neptune, respectively. For Neptune, 
a large fraction of the planet was visible through the slit except 
for the E and W limbs. 

The employed grism (EFOSC2 grism#5) covers the wave- 
length range of 5200 - 9350 A and provides a spectral resolution 
of 12.8 A for a 1" wide slit. The data were recorded with a 2kx2k 
CCD (ESO CCD#40) with a spatial scale of 0.157" pixel" 1 and 
a wavelength scale of 2.06 A per pixel. 

The effective seeing for our spectropolarimetric data was 
about 1" as derived from the width of the spectra of the stan- 
dard stars. 

Our polarimetric data are expressed in form of the normal- 
ized Stokes parameters Q/I and U/I. There is Q/I = (7|| - 
I±W\\ + I±) = (7o - hoWo + bo), where /y = I and I ± = 7 90 
are the intensities with polarization direction parallel and per- 
pendicular to the slit, respectively. The Stokes parameter U/I — 
(Z45 — /i35)/(/45 +/135) describes polarized light in the 45°/135° di- 
rection, where positive angles are rotated counterclockwise. The 
polarization position angles are defined relative to the slit orien- 
tation. 

The linear polarization was measured in a standard way (e.g. 
Tinbergen & Rutten 1992) with sets of four exposures taken with 
half-wave plate position angles of 0°, 22.5°, 45° and 67.5°, re- 
spectively. Three sets were taken for each slit orientation and 
both planets in order to enhance the signal-to-noise. One sin- 
gle exposure was 30 sec for Uranus and 50 sec for Neptune. 
Polarimetric standard stars were observed with the same instru- 
mental setup. Exposures of a helium-argon lamp provided the 
wavelength calibration. 

Exposures with the half- wave plate angles at 0° and 45° yield 
Stokes Q and the other two yield Stokes U. One Q/I spectrum 
results from the combination of 4 planet spectra according to 



Q 
I 



R-l 
R + l 



with R z 



/||(0 o )/^(0°) 
/|l(45 )// i (45°) 
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and equivalent for U/I with the I\\ and I ± taken at half-wave plate 
angles of 22.5° and 67.5°, respectively. 

The polarization reduction was performed with the midas 
software package similar to the description in paper 1. Important 
in the data reduction procedure is that the different planetary 
spectra, e.g. 7||(0°), 7 ± (0°), /||(45°) and 7 ± (45°), are centered in 



the same way and with an accuracy of about 1/10 of a pixel 
in spatial direction. This high precision is required because a 
small misalignment would produce a spurious polarization sig- 
nal (say positive) at one limb and an opposite polarization sig- 
nal (negative) at the other limb of the long-slit spectrum. For 
an accurate alignment it has to be taken into account that the 
Wollaston has some dispersive power. The resulting ordinary and 
extra-ordinary spectra are therefore slightly bended with a cur- 
vature which differs between the two beams. The alignment of 
the long-slit spectra of the planets was achieved with the help 
of the standard star observations from which the wavelength de- 
pendent spatial position of the spectra could be determined with 
sufficient precision. 

The polarization zero point of the instrument was de- 
rived and corrected from an unpolarized (HD 14069) and 
a polarized (BD +25°727) standard star. The instrumental 
polarization was found to be less than 0.2 %. According 
to the EFOSC2 documentation the instrument polarization 
should be less than 0.1 % for the central parts of the CCD 
(http://www.ls.eso.org/lasilla/sciops/3p6/efosc/). The polariza- 
tion angle calibration should be accurate to about 9 « ±2°. 

The CCD introduces for A > 7000 A an interference pattern 
with a full amplitude of about 5 % at 7500 A and about 10 % at 
8200 A of the mean intensity level. The differential data reduc- 
tion for the polarization reduces this fringes to a full amplitude 
level for Q/I of Ap * 0.3 % at 7500 A and * 0.5 % at 8200 A. 
In the intensity and polarization spectra the fringes are visible as 
quasi-periodic pattern with a periodicity of about 15 - 20 A. We 
have tried to reduce the fringe pattern with the available flat field 
calibrations. Some improvement was achieved for the intensity 
spectrum but not for the polarization spectrum. Spectral binning 
turned out to be much more efficient. This averages out the peri- 
odic pattern and produced much smoother intensity and polariza- 
tion spectra. As the spectral features in Uranus and Neptune are 
very broad the loss in spectral resolution was acceptable. With 
this method the spurious fringe pattern vanished almost com- 
pletely and in addition a significantly improved S/N per bin was 
obtained for the polarization spectrum. Even with no fringes a 
similar spectral binning would have been necessary for an analy- 
sis of the polarization spectrum in order to have sufficient signal- 
to-noise. For these reasons we have binned all spectropolarimet- 
ric data into 30 A bins. 



3. Spectropolarimetric structure of Uranus 

For Uranus we have taken 5300 to 9300 A long-slit spectropo- 
larimetry with the slit oriented in N-S and E-W direction in the 
celestial coordinate system (Fig. Q]). The N-S slit position ex- 
tends roughly along the equator of Uranus. Henceforth, measure- 
ments taken with this slit position are called "equatorial" with 
the definition of West and East limbs or longitudes on the planet 
as indicated in Fig.Q] 

The E-W slit extends from the South pole over the equator to 
northern latitudes roughly along the central meridian of Uranus. 
Measurements from this slit orientation will be referred to as 
"meridional". 

Figure Q] illustrates the exact slit positions and the selected 
regions: four limb regions East, West, South and NortrQ (0.628" 
long black boxes), two central slit regions (0.942" long white 



1 From here on East, West, South and North refer always to the plan- 
etary coordinates. 
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Fig. 1. Uranus i-band images with the slit in celestial N-S (left) 
and E-W (right) orientation. N is up and E is left. The black 
boxes indicate the "limb" sections of the long-slit spectropo- 
larimetry and the white boxes the "central" regions. Also shown 
are the positions of the South pole and the equator with the West 
limb indicated according to the planetary coordinates defined by 
Seidelmann et al. ( 120021 1. At the time of our observations the 
diameter was 3.51", the phase angle 2.8°, the position angle of 
the South pole was 78° and the distance of the South pole from 
the disc center was 1 .62". (The planet parameters are taken from 
"The Astronomical Almanac 2003".) 



boxes), as well as the two spectropolarimetric signals obtained 
by integrating over 7.85" along the equatorial and the meridional 
slit. These selected regions are denoted as limb, central part and 
total slit. 

Already our ground based images reveal a strong asymme- 
try between the brighter South pole and darker northern latitudes 
which are well resolved as band structures in HST images (e.g. 
Karkoschka 1200 1 1 Rages et al. [2004), but no strong longitudi- 
nal features are seen for Uranus (see also press release STScI- 
PRC2004-05). 

3. 1 . Spectroscopic structure and spectral albedo A(A) 

Spectroscopically (Fig. [2] top panel), our data show mainly the 
strong methane (CH4) absorption bands, which are well known 
for Uranus (e.g. Baines et al. 1 1 9861 Karkoschka 1994). Also 
clearly visible is the telluric absorption due to the A-band at 
7590 A. Further, the enhanced limb brightening in the methane 
bands is well visible in our data and illustrated (in Fig.|2]i by the 
normalized limb to albedo ratio. 

An albedo spectrum is constructed from our data. For this we 
averaged first the total slit spectra in equatorial and meridional 
direction and then normalized the result to the geometric albedo 
spectrum of Karkoschka ( 1994) multiplying our spectrum with a 
smooth correction function which accounts for the instrumental 
efficiency curve. This procedure yields a well defined full disk 
albedo spectrum A (A) as plotted in Fig. [2] 

Spectral regions were selected in order to quantify the albedo 
and polarization. The selected regions, which represent strong 
absorption bands (Ml - M7) and the higher albedo regions (CI 
- C5) in between, are indicated in Fig.|2]and listed in TableQ] 

We found a clear difference of the spectral structure of the 
limb when compared to the albedo spectrum A (A). This is shown 
in Fig. [2] where the normalized ratio between the average of the 
West and East limb spectra and the albedo spectrum is plotted. 
The ratio is enhanced in all CH4 bands indicating that absorption 
is less deep for the limb regions. This is equivalent to the state- 
ment that deep absorption bands of Uranus show a limb bright- 
ening effect as observed e.g. by HST (see Karkoschka 2001). 




6000 7000 8000 9000 

wavelength 



Fig. 2. Spectropolarimetry of Uranus for the equatorial slit. Top 
panel: geometric albedo spectrum A(A) with identifications of 
selected spectral bands. Second panel: normalized ratio between 
the limb intensity spectrum and the geometric albedo spectrum 
A(A). Third panel: normalized Stokes Q/I polarization spectra 
for the average of the western and eastern limb regions (solid 
line), the central part (dotted line), and the total slit (dashed line). 
Bottom panel: same as third panel but for Stokes U/I. 



3.2. Spectropolarimetry for the equatorial slit 

Spectropolarimetry for the equatorial slit may be considered as 
representative for the "average" planetary disk of Uranus be- 
cause the bright South pole region and the faint northern lati- 
tudes are not included in these observations. The peculiarities of 
these special regions are discussed in the following subsection. 

The polarimetric signal of the equatorial cut through Uranus 
is plotted in Figure [2] as normalized Stokes Q/I and U/I spectra 
(Pq(A) and pu(A) respectively). There exists in our data no sig- 
nificant difference between the polarization signal of the western 
and eastern limbs and therefore only the average of the two limb 
polarization spectra is plotted. Fig.|2]includes further the polar- 
ization spectra of the central part and the average for the entire 
slit. 

The data show, that the 2/^-polarization is positive and high 
at the limbs and close to zero for the central region. Essentially 
no signal is visible for the £///-data, indicating that the orienta- 
tion of the polarization is everywhere parallel to the slit equiv- 
alent to a limb polarization perpendicular to the limb. The U/I- 
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spectra give a good representation of the noise errors in the data, 
which are particularly large in the deep absorption bands at the 
longest wavelength A > 8500 A, where the photon statistics are 
poor. 

The ^//-polarization spectrum for the equatorial limbs dis- 
plays a clear overall decrease of the polarization level with wave- 
length. Further, the limb polarization is enhanced in the deep 
CH4 absorption bands when compared to the adjacent contin- 
uum or spectral regions with higher albedo. 

TableQ]gives the measured polarization for all four limbs for 
selected wavelength regions. The listed values are flux weighted 
average values (Q/I) = SQ/S7, where "EI and EQ are summed 
over intensity flux and polarization flux respectively for the se- 
lected spatial bin and wavelength interval. 

In the methane bands the limb polarization decreases from 
about 2.6 % in the band at 5400 A to about 0.9 % in the broad, 
deep band centered at 8900 A. The measured limb polariza- 
tion for the continuum (or high albedo regions) is about 2 % at 
5500 A, and less than 0.4% for the C5-region at 8270 A, except 
for the northern limb where the polarization for C5 is signifi- 
cantly higher. 

It should be noted, that the level of polarization in Fig.|2]and 
Table [TJ depends on the seeing conditions and the size and posi- 
tion of the spatial bin, selected for the averaging. However, the 
relative wavelength dependence of polarization for a given limb 
section depends very little on spatial position, and therefore on 
the selected slit region and the seeing conditions. Uncertainties 
in the polarization values of Table Q] due to photon noise are 
about A{Q/I) = +0.10%, except for the narrow interval M3 
and the longest CH4 wavelength intervals M6 and M7, where 
A(Q/I) = +0.20%. 



3.3. Peculiarities for the meridional slit 

Uranus, close to its equinox, which is on the 7 th of December 
2007, shows currently an intensity asymmetry from the bright 
South pole to the darker northern latitudes which is particularly 
strong in the methane bands. This intensity asymmetry is also 
obvious in the HST-images of Karkoschka ( 120011 1 or Rages et 
al. (2004) and in Keck observations obtained with an adaptive 
optic system (Sromovsky & Fry |2005l l. 

Also in the spectropolarimetric signal we found some clear 
peculiarity for the northern latitudes. Figure[3]compares the Q/I 
limb spectropolarimetry of the northern latitudes with the aver- 
age of the three other limb regions, East, South and West which 
are very similar and can be taken together. The northern latitudes 
show for wavelength longer than about 6000 A no polarization 
minimum for spectral regions with high albedo. Thus the nor- 
malized polarization spectrum Q/I for the northern latitude limb 
is essentially featureless and displays just a steady decrease in 
polarization with wavelength (see also Table[3]l. 

According to our data, the darker northern latitudes show, 
when compared to the other limbs, for A < 6000 A a higher over- 
all polarization and for A > 6000 A a higher polarization for the 
inter band regions. We have carefully investigated whether this 
result could be due to a misalignment or another spurious effect. 
We found that the maximum error in misalignment for the po- 
larimetric data reduction produces an effect on the order of only 
0.05. Based on such tests, we conclude that the described Stokes 
Q/I polarization peculiarities of the northern limb are real. 



Table 1. Limb polarization (Q/I) for Uranus for the West, East, 
South and North limb and selected wavelength intervals (column 
2). Column 3 gives the averaged albedo A(A) for the correspond- 
ing interval. 
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1.88 
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1 Q1 
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0.41 


1.40 
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1.85 


C4 


7442 - 7565 


0.32 


0.89 


0.88 


0.75 


1.47 


C5 


8215 - 8320 


0.23 


0.27 


0.33 


0.20 


1.10 


Ml 


5368 - 5459 


0.37 


2.68 


2.63 


2.20 


3.10 


M2 


5732 - 5788 


0.40 


2.41 


2.26 


1.96 


2.60 


M3 


6180 - 6200 


0.13 


2.45 


2.52 


1.96 


2.29 


M4 


6627 - 6699 


0.21 


1.84 


1.82 


1.51 


1.93 


M5 


7174-7358 


0.06 


1.55 


1.70 


1.20 


1.40 


M6 


7933 - 8034 


0.05 


1.27 


1.25 


0.96 


1.37 


M7 


8369 - 8453 


0.06 


0.88 


1.03 


0.66 


1.21 
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Fig. 3. Comparison between Q/ /-polarization spectra of the 
northern limb of Uranus (thick line) and the average of the other 
three limbs (thin line). 

3.4. Fitting the limb polarization spectrum 

We found that the Q/I = pq(A) polarization spectra for Uranus 
show a clear overall decrease of the polarization level with wave- 
length. Further, the polarization is enhanced in the deep CH4 ab- 
sorptions bands, when compared to the adjacent higher albedo 
spectral regions, except for the red part of the northern limb 
spectropolarimetry. 

In order to quantify the wavelength and albedo dependence 
of the limb polarization, we fit the polarization spectra with a 
simple linear relation with respect to the wavelength A and the 
albedo A(A): 

p Q (A)[%] = Cl + c 2 ■ A[pm] + c 3 ■ A(A) . (2) 

The parameters c\, c 2 , and C3 can then be determined with a 
least square procedure. Particularly good fits were obtained for 
the spectral region from 5300 A to 7500 A (see Fig . |4j» . 

The best fit to the equatorial limb polarization for the wave- 
length region 5300 - 7500 A is: 

p Q (A)[%] = 8.4 - 9.1 ■ A[pm] - 2.2 ■ A(A). (3) 
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Fig. 4. Fits to the equatorial limb polarization spectrum com- 
pared to the measured spectrum (thin line). The thick line is a 
best fit for the spectral region 5300 - 7500 A and the dotted line 
for 7700 -9300 A. 



The standard deviation of the fit is less than 0.1%. The 
largest deviations are due to the telluric absorption bands at 
6870 A (B-band) and 7200-7300 A (water absorption) which 
were not corrected in the albedo spectrum. Due to the uncor- 
rected A-band absorption at 7590 A, the upper boundary of the 
least square fitting procedure was fixed to 7500 A. 

The fit to the longer wavelength region 7700 A to 9300 A 
yields a flatter wavelength dependence and a stronger albedo de- 
pendence (ci = 6.0, C2 = -5.6 and C3 = -4.5). 

3.4.1 . Spatial sampling and full slit signal 

According to a simple analytic model for Rayleigh scattering 
atmospheres (see e.g. paper 1; van de Hulst 1980) the limb polar- 
ization feature is expected to be a narrow, strongly peaked fea- 
ture along the limb. In our data this narrow peak is not resolved 
due to the seeing limited (« 1 arcsec) resolution. Therefore, 
the apparent strength of the limb polarization pq(A) depends 
strongly on the not so well defined seeing conditions during the 
observations. In addition we have to consider also the slit width 
which defines the spatial region over which the intensity and 
Stokes fluxes are sampled. For Uranus we employed a narrow 
slit (0.5 arcsec) which caused only a small degradation of the 
polarization signal. 

Choosing a radial bin for the quantification of the limb polar- 
ization is an arbitrary procedure. For this reason it is reasonable 
to employ the polarization signal integrated over the entire slit 
to quantify the limb polarization. A fit according to Eq. (O to 
the total slit spectropolarimetry of the equatorial slit as given in 
Fig-Elyields the parameters c\ = 3.6, C2 = -3.8, and C3 = -1.3 
(again for the spectral region 5300 - 7500 A). 

In Sect. 14.2.11 we derive a correction factor which accounts 
for the seeing and the spatial sampling of our observations. 



4. Disk profiles for Uranus 

The spectropolarimetric data can be employed to investigate the 
intensity and polarization profiles through the disk of Uranus for 
any wavelength interval covered by our observations, a sort of 
center to limb profiles. Since the slit is much longer than the 
planet's diameter our data provide profiles roughly along the 
equator (equatorial slit) and profiles from the South pole to high 
northern latitudes (meridional slit). 



Again, we focus our analysis on the spectral regions CI - C5 
(continuum /high albedo regions) and Ml - M7 (strong methane 
bands) as defined in Table Q] 

Figures [5] and [6] show examples for disk profiles for the Cl- 
region around 5850 A and the M2 methane band at 5760 A. The 
flux of the disk profiles was normalized at the disk center Iq = 
I(x = 0") = 1 and also the Stokes flux Q is expressed relative to 
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Fig. 5. Uranus equatorial (solid) and meridional (dashed) disk 
profiles for the continuum (high albedo) wavelength range 
5827 A to 5883 A (CI). The panels give the flux /(jc)// (top), 
the normalized Stokes parameter Q{x)/I{x) (middle) and the 
Stokes flux Q(x)/I Q (bottom). 



4. 1 . Disk profiles for the intensity 

The intensity profiles show some pronounced differences be- 
tween the equatorial and the meridional directions as well as 
between different spectral regions. All equatorial profiles are es- 
sentially symmetric while the meridional profiles are more or 
less asymmetric, showing a brighter South pole when compared 
to the northern latitude (see Fig. |6). This asymmetry is partic- 
ularly strong in the deep methane bands. In this section we fur- 
ther present a relationship between intensity profile structure and 
albedo for both slit orientations. 
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Fig. 6. Uranus equatorial (solid) and meridional (dashed) disk 
profiles for the methane band centered on 5760 A (M2). The 
panels give the flux I(x)/Iq (top), the normalized Stokes param- 
eter Q{x)/I(x) (middle) and the Stokes flux Q{x)/Iq (bottom). 



The equatorial cuts have a Gaussian-like profile for the con- 
tinuum / high albedo bands (Fig. [5]) and flat top profiles for the 
deep absorptions (Figs.|6]and|7]). The flat top profiles can be ex- 
plained as not well resolved limb brightening effect. 

These intensity features of Uranus are also clearly present 
in imaging studies based on high resolution HST observations 
(e.g. Karkoschka l200Tl Rages et al. l2004j l. Compared to the HST 
data, we can only see the strongest spatial features due to the 
much lower spatial resolution of our data. However, the spec- 
tropolarimetric data yield disk profiles for any spectral region in 
the covered wavelength range and therefore our data are partic- 
ularly well suited to investigate the systematic behavior of the 
asymmetry and the limb brightening for Uranus. 

In order to quantify the asymmetry between South and North 
and the limb brightening effect we construct first a symmetric 
disk profile as reference profile 7 re f (x) for Uranus by averaging 
all 5 equatorial continuum profiles (CI - C5) and their mirrored 
profiles (mirrored at x = 0")- As next step, we determine the 
relative deviation 6(x) for the equatorial and meridional profiles 
from the reference profile by calculating 



S(x) = 



I(x) - 7 ref (x) 



(4) 



Table 2. Uranus meridional asymmetry parameter A mer and limb 
brightening parameter Y, eq at the equatorial regions for spectral 
features CI - C5 and Ml - M7 as defined in Table Q] 



feature 






CI 


-0.017 


-0.082 


L2 


-0.030 


-0.019 


C3 


-0.008 


0.019 


C4 


0.038 


0.048 


C5 


0.102 


0.035 


Ml 


0.031 


0.070 


M2 


0.079 


0.039 


M3 


0.095 


0.375 


M4 


0.141 


0.149 


M5 


0.160 


0.529 


M6 


0.102 


0.467 


M7 


0.156 


0.341 



Due to our normalization, the deviation in the middle of the 
profiles is for all our profiles 5(x = 0") = 0. Strong signals in 
6(x) are obtained for x » ±1". We calculated the averages of the 
deviations (6{x)) for all four planetary orientations on the disk 
selecting profile regions around x ~ 1". The selected regions 
correspond to the "grey" boxes in Fig.Q] defined by: 



5e,n = (8{x)) 



for -1.41" < x < -0.47' 
for +0.47" <x< +1.41' 



Based on this we calculated the meridional asymmetry pa- 
rameter: 



turner - <$s - 8\ 



and the equatorial limb brightening parameter: 



(5) 



(6) 



which are given in Table [2]for the wavelength intervals CI - C5 
and Ml - M7. Large values indicate a relatively strong asymme- 
try between South and North or a relatively strong limb brighten- 
ing. Figure Q illustrates the asymmetry for the central meridian 
and the equatorial limb brightening for the most extreme case, 
the M5 absorption band at 7270 A. 

It should be noted that the absolute values of A mer and T, eq 
contain little diagnostic information, because they depend on the 
seeing and the size of the chosen spatial interval. 

More interesting information is obtained from the relative 
strength of the A mer and T, e „ parameters for different wavelength 
bands. They are all affected in a similar way by the limited spa- 
tial resolution. Our data show a tight anti-correlation between 
the meridional asymmetry parameter A mer and albedo A as well 
as between limb brightening T, eq and A which is illustrated in 
Fig.! 

The HST observations of Uranus described by 
Karkoschka (2001) also clearly show that the intensity 
asymmetry and limb brightening are particularly strong in the 
methane bands. However, our long-slit spectroscopy allows 
to select any particular wavelength band for a center-to-limb 
profile analysis. To our knowledge, the dependence of the 
asymmetry and limb brightening effect on the albedo has not 
been described previously. The asymmetry between South pole 
and northern latitudes is of course a seasonal effect of Uranus. 
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Fig. 7. Disk intensity profile for the methane absorption M5, 
showing the strongest asymmetry between South and North (dot- 
ted line) and equatorial limb brightening (solid) effect. For com- 
parison the symmetric reference profile 7 le f is plotted in addition 
(Gaussian-like). 
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Fig. 8. Uranus meridional asymmetry parameter A mer vs. albedo 
(top) and limb brightening 1 eq vs. albedo (bottom) for high 
albedo spectral regions CI - C5 (□) and methane absorption 
bands M1-M7 (o). 



To investigate the properties of the scattering layers in 
Uranus and Neptune via the limb brightening effect, it would be 
desirable to have long-slit spectroscopy with a spatial resolution 
as offered by HST. 

4.2. Disk profiles for the polarization 

Disk profiles for the polarization are obtained for the normal- 
ized Stokes parameter Q(x)/I(x) (short Q/I) and the Stokes flux 
Q(x)/Iq (short Q), where positive values indicate a polarization 



parallel to the slit or perpendicular to the limbs. A few examples 
are plotted in Figs. [5] and [6] 

As already described in the previous sections, the polariza- 
tion is low in the center of the disk and high at the limb at all 
wavelengths of our spectrum. Typically Q/I increases steadily 
with radius r = \x\, and reaches a constant value around r ~ 2", 
where the photon noise starts to dominate the measurements. 
The Stokes flux Q increases from r = 0" until about r » +1.5", 
and decreases together with the intensity profile further out to 
zero. This behavior is in agreement with the imaging observa- 
tions described in paper i. 

For the equatorial slit observations all Q/I and Q disk pro- 
files are essentially symmetric, and they have roughly the same 
spatial structure (Figs. [5] and |5). Only the strength of the polar- 
ization signal shows a clear wavelength dependence. 

For the meridional slit there is an obvious asymmetry in the 
Q/I polarization signal. For most spectral regions the Q/I polar- 
ization is typically 30 % higher at the northern limb than at the 
South limb (Tab.Q] Figs. [5] and |6}. The asymmetry is particularly 
strong for the inter band regions C4 and C5. This effect is also 
described in Sect. 13.31 where the spectropolarimetric signal for 
the southern and northern limb is compared. 

Surprisingly, the Q-fiux profiles in meridional direction are 
rather symmetric (Figs.[5]and[6]l. Thus, the higher intensity I at 
the South pole, combined with the lower normalized polarization 
Q/I, yield together a Stokes flux Q at a similar level as for the 
darker, but higher polarized northern latitudes. Whether this is a 
fortuitous coincidence or an indication of a planet-wide homo- 
geneous haze layer producing everywhere more or less the same 
limb polarization flux remains to be investigated. 

4.2.1. Seeing corrected limb polarization parameters for 
Uranus 

For comparison with models and future observations of Uranus, 
it appears useful to derive well defined parameters for the 
strength of the limb polarization. In particular, such values 
should be corrected for the polarization cancellation introduced 
by the seeing, and the averaging effects due to the finite slit 
width. 

For this reason we first determine the flux weighted mean 
polarization for a line (slit) through the center of the planet 
{Q/I)\me- This is calculated according to (Q/I)n ne = "LQ/YA, 
where "LQ and 1,1 are the polarization flux Q{x) and flux I(x), 
summed up along the entire slit (identical to the procedure in 
paper i). The slit integrated polarization {QII)\\ ne is of course 
significantly lower than the values derived for a small limb re- 
gion given in Table[T]or shown in Fig.[2]for the equatorial slit. 

The seeing introduced by the Earth's atmosphere causes a 
polarization cancellation effect. As explained in detail in paper 
i, the opposite sky polarization components, e.g. the +2sky com- 
ponents at the equatorial limbs, and the -Qsky components at the 
northern and southern limbs, overlap and cancel each other if the 
seeing is substantial, say seeing > 0.5 Rpianet- A second effect of 
the seeing is that the signal from regions adjacent to the slit spill 
into the slit. This averages down the measurable limb polariza- 
tion. In addition, there is also a reduction of the integrated limb 
polarization due to the finite slit width. For Uranus this is a small 
effect, because the data were obtained with a very narrow slit of 
0.5". 

For our observations the seeing-limited resolution was de- 
termined to be 1.09" at 6000 A and 1.00" at 8000 A taken from 
observations of standard stars observed with the same setup. 



8 



F. Joos and H.M. Schmid: Limb polarization of Uranus and Neptune. II 



However, the seeing conditions can change on short time scales, 
introducing uncertainties on the order of +0.1" to +0.2". 

This effect of polarization cancellation due to the limited res- 
olution was modeled in paper i, and seeing correction factors 
for the imaging polarimetry were determined. We have repeated 
these calculations for observations through a 0.5" wide slit and 
obtained essentially the same seeing dependence as for the imag- 
ing polarimetry (see Fig. 1 1 in paper i). From these simulations 
we found that the measurable polarization is reduced by a factor 
0.96 due to the width (0.5") of the slit and a factor of 0.76 due to 
the seeing of 1" when compared to the intrinsic limb polariza- 
tion integrated along a slit through the disk center. 

Thus, in Table [3] we give corrected values for the intrinsic 
limb polarization (Q/I)^ for Uranus by multiplying the mea- 
sured values (Q/I)n ne with a factor of 1.37. The uncertainty in 
this correction factor is about +0.15 mainly due to the not pre- 
cisely defined seeing. For the intrinsic limb polarization other 
error sources like photon noise or the calibration of the polariza- 
tion data are less important. 

Table 3. Line (slit) integrated, seeing corrected polarization 
(Q/I)\°™ for the equatorial and meridional slits through the cen- 
ter of Uranus for the selected spectral features (Table Q]). The 
listed values are corrected (measured values xl.37; see text) for 
the seeing and the employed slit width of our observations. 



feature (QIDnZ [%1 

equatorial South North 



CI 


0.90 


1.02 


1.19 


C2 


0.80 


0.99 


1.13 


C3 


0.66 


0.90 


1.05 


C4 


0.41 


0.57 


0.89 


C5 


0.18 


0.24 


0.67 


Ml 


1.47 


1.46 


1.85 


M2 


1.26 


1.40 


1.55 


M3 


1.61 


1.46 


1.57 


M4 


1.01 


1.13 


1.28 


M5 


1.15 


0.95 


1.04 


M6 


0.88 


0.75 


1.03 


M7 


0.71 


0.46 


1.08 



The same correction factor can also be applied to the fit to the 
equatorial slit-integrated spectropolarimetric signal determined 
in Sect. 13.4.11 Thus, the corrected fit for the intrinsic equatorial 
full slit polarization for the spectral region 5300 to 7500 A can 
be given by: 

/? e (^) = 4.9-5.2^[ / um]-1.8A(T). (7) 

5. Spectropolarimetric structure of Neptune 

Spectropolarimetry of Neptune was obtained with a North- 
South oriented slit and almost the same instrument set-up as for 
Uranus. The only difference was the slit width of 1.5", which 
was much broader than for Uranus. We also define three slit 
regions for Neptune, a North limb and South limb section of 
0.625" each, a central region of 0.942" and an integrated "total 
slit" of 4.71" length (Fig. |9]). 

At the time of our observations Neptune had a diameter of 
2.24". Compared to this size the width of the slit of 1.5" and 
the seeing conditions of about 1" are large and not well suited to 
resolve spatial features on the planetary disk. Almost the entire 



planet was contained in the slit except for the eastern and western 
limbs. With higher spatial resolution observations it was shown 
that Neptune displays also a limb brightening effect similar to 
Uranus (e.g. Hammel et al. 1 19891 Baines & Hammel [1994l ). 




Fig. 9. Intensity image of Neptune in the i-band with North- 
South slit position. North is up and East is left. The South pole 
is marked with a dash. The grey stripe indicates the position of 
the slit. The black boxes show the "limb" area and the white box 
the "central part". At the time of our observations the diameter 
was 2.24", the phase angle 1.7°, the position angle of the South 
pole 169° and the distance of the South pole from the center of 
the disc was 0.97". 

Despite the bad resolution of our Neptune observations a 
clear polarization signal is detected. This is shown in Fig. [10] 
which presents the spectropolarimetry of Neptune in a similar 
form as for Uranus (Pig. |2j. The albedo spectrum A(A) was con- 
structed with a normalization to the Neptune albedo spectrum of 
Karkoschka (1994). Spectroscopically, Neptune is very similar 
to Uranus (see Karkoschka l 19981 1. 

The spectropolarimetric signal of Neptune for the average of 
the northern and southern limbs is qualitatively also very sim- 
ilar to Uranus (Fig ITOb. The Q/I polarization spectrum is high 
and positive at the limbs, or parallel to the slit, but is low in the 
disk center, and essentially zero for the f/-direction. Q/I shows 
an overall decrease in polarization towards longer wavelength 
and an enhanced polarization in the strong methane bands, when 
compared to the adjacent continuum or higher albedo regions. 
The level of the measured limb polarization in Neptune is sig- 
nificantly lower than for Uranus because the planetary disk is 
less well spatially resolved. 

The spectropolarimetry for the central region shows a weak 
and negative Q/I polarization spectrum. This can be easily under- 
stood as polarization signal from the eastern and western limb 
which spills into the central slit region. As the eastern and west- 
ern limb polarizations are oriented perpendicular to the limb, 
they produce a perpendicular (negative) signal for the North- 
South slit. The slit integrated polarization is small because of 
the small spatial resolution. For spatially unresolved observa- 
tions with a slit which is wider than Neptune, a total net polar- 
ization close to zero is expected, because the centro-symmetric 
polarization structure of Neptune (see paper i) would produce an 
almost perfect cancellation of positive and negative polarization 
features. 

Parameters for the limb polarization of Neptune are given 
in Table |4] Because Neptune is spectroscopically very similar 
to Uranus, the same spectral features were selected. The mean 
albedo for the spectral intervals was determined from the albedo 
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Fig. 10. Spectropolarimetry of Neptune for the slit in North- 
South direction. Top panel: geometric albedo spectrum A(A). 
Middle panel: normalized Stokes Q/I polarization spectra for 
the average of the North and South limb regions (solid line), the 
central part (dotted line) and the entire slit (dashed line). Bottom 
panel: Same as the middle panel but for Stokes U/I, 

spectrum A(A). The limb polarization is only given for the av- 
erage of the northern and southern limb regions. The two mea- 
surements are averaged because no significant difference is de- 
tectable between the North and the South limb. 

Because the spatial resolution of our observation is poor, it is 
not possible to deduce reliable parameters for the intrinsic limb 
polarization of Neptune. However, the imaging polarimetry pre- 
sented in paper i indicates that the intrinsic limb polarization of 
Neptune is about 1.5 times higher than that of Uranus. 



Table 4. Limb polarization (Q/I) (in [%]) for Neptune for the 
average of the North and South limbs, (Q/I)ns an d f° r the total 
meridional slit (Q/I)t, for selected wavelength intervals. Column 
2 gives the averaged albedo A(A) for the corresponding interval. 



feature 


albedo 


(Ol 'Dni 


(QIDt 


Cl 


0.49 


0.85 


0.37 


CO 




u.ou 


U.JU 


C3 


0.30 


0.77 


0.34 


C4 


0.22 


0.62 


0.29 


C5 


0.14 


0.38 


0.23 


Ml 


0.32 


1.22 


0.52 


M2 


0.33 


1.09 


0.48 


M3 


0.11 


1.08 


0.47 


M4 


0.16 


0.97 


0.42 


M5 


0.06 


0.83 


0.38 


M6 


0.05 


0.67 


0.33 


M7 


0.05 


0.67 


0.33 



This indicates that the wavelength and albedo dependence of the 
limb polarization is very similar for Neptune and Uranus. 

For the red end of the North-South spectropolarimetry (7700 
- 9300 A) the derived parameters are c\ = 1.75, C2 = —1.15 and 
C3 = -2.75. For the longer wavelength region the wavelength 
dependence is flatter and the albedo dependence stronger. 




6000 



7000 8000 
wavelength 



9000 



Fig. 11. Fits to the North-South limb polarization spectrum of 
Neptune compared to the measured spectrum (thin line). The 
thick line is a fit for the spectral region 5300 - 7500 A, and the 
dotted line for 7700 - 9300 A. 



5.1. Fitting the polarization spectrum of Neptune 

We present also a fit to the limb polarization spectrum of 
Neptune for a quantitative description. The same procedure is 
applied as described in Sect. l3.4l This yields for Neptune for the 
average of the North and South limb polarization spectrum and 
the wavelength range of 5300 - 7500 A the fit: 

p Q (A)[%] = 3.35 - 3.35 ■ A[jxm] - 1.10 -A(A). (8) 

It is interesting to note that the fit parameters for Neptune be- 
have like ci m —ci and c\ » -3c3. These are almost identical 
proportionalities as derived for the limb polarization in Uranus. 



6. Discussion and Conclusions 

6.1. Comparison with the imaging polarimetry ofpaperi 

The R, i and z-band imaging polarimetry of Uranus and Neptune 
reported in paper i reveals the same overall polarization behav- 
ior as the spectropolarimetric data presented in this work. Both, 
imaging polarimetry and long-slit spectropolarimetry, show for 
Uranus and Neptune essentially no polarization in the disk center 
and a strong limb polarization with a position angle perpendicu- 
lar to the limb. 
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Already the polarimetric images taken in a few, mainly 
broad-band, filters define a general trend, that the limb polariza- 
tion is lower for longer wavelengths. However, the wavelength 
dependence of the limb polarization is described in much more 
detail with the long-slit spectropolarimetry presented in this pa- 
per. From spectropolarimetry it is obvious that there is not only 
a decrease in the limb polarization with wavelength but also a 
very tight anti-correlation with the albedo. 

Our spectropolarimetry and the imaging polarimetry of 
Uranus and Neptune yield further results which are complemen- 
tary. 

Very interesting is the spectropolarimetric result that the rel- 
atively dark northern latitudes have a higher level of polarization 
Q/I when compared to the other limbs. With the imaging po- 
larimetry it was not possible to find this result due to the delicate 
alignment requirements for the images with opposite polariza- 
tion in the data reduction process. 

The imaging polarimetry was taken under better seeing con- 
ditions (0.8") and the spatial resolution was not further degraded 
by a wide slit, as it was for Neptune in our spectropolarimetry. 
Therefore, the imaging polarimetry allowed to derive a seeing 
corrected limb polarization average for both planets, indicating 
that the intrinsic limb polarization for Neptune is about a factor 
of S3 1.5 higher than for Uranus. 

A good quality test of the data analysis provides the com- 
parison of the derived intrinsic limb polarization of Uranus from 
imaging polarimetry and spectropolarimetry. For this we have to 
convert the derived values for the disk averaged limb polariza- 
tion (Q r /I) from imaging polarimetry into the slit (line) aver- 
aged limb polarization (Q/I)u n e as obtained from the spectropo- 
larimetric data. This conversion depends somewhat on the radial 
dependence of the intensity and the limb polarization. However, 
the simple model results presented in paper i indicate that a con- 
version relation (Q r /I) ~ 1-6 (Q/I)\i ne is adequate for the ex- 
pected center to limb profiles. 

Thus, the converted parameters for the intrinsic limb po- 
larization of Uranus from the filter polarization is (Q/I)™™ = 
0.53 % for the %73-filter (albedo=0.32), 0.59 % for the % 29 - 
filter, corresponding to M5, (albedo=0.05), 0.46 % for the i-filter 
(albedo=0.19) and 0.27 % for the z-filter (albedo=0.12). 

In Fig.[T2lthe limb polarization (Q/I)^ and the albedo from 
the filter polarimetry are compared with the spectropolarimetry. 
The broad band filters i and z lie between the values for the red- 
dest wavelength intervals C5 and M6 and M7 as expected for a 
broad band average. Also for the narrow band filter ^73 located 
in wavelength between M4 and C3 the agreement is quite good. 
Only for the A^9 filter the measured polarization from the filter 
polarimetry is significantly lower than the spectropolarimetric 
results (M5). As described in paper 1, the accuracy of the mea- 
surement in this filter is particularly low due to the low signal- 
to-noise. Thus, it may be concluded that the agreement between 
imaging polarimetry and spectropolarimetry is good, except for 
the low quality /1729-firter observations. 

This gives confidence that the analysis and the applied seeing 
corrections are consistent for the two data sets. 

6.2. Limb polarization and atmosphere models 

The properties of the limb polarization for Uranus and Neptune 
can be summarized as follows: higher polarization for shorter 
wavelength, enhanced polarization in the methane bands when 
compared to the continuum or inter band regions, and in Uranus 
a higher continuum/inter band polarization for the darker north- 
ern latitudes when compared to the bright South Pole. 
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Fig. 12. Seeing corrected, full slit polarization for the equato- 
rial slit of Uranus vs. albedo from Table [3] The circles are the 
high albedo regions (CI - C5) and the triangles are the methane 
absorption bands (Ml - M7). The dotted line gives the full slit 
polarization at 6000 A, and the dashed line, at 7000 A, according 
to Eq. (0. The filled squares denote the measured polarization 
with the four filters from paper 1 (broad band filters i and z and 
the narrow band filters centered at 673 and at 729 nm). 



An enhanced polarization in the strong methane bands is also 
observed for Jupiter (Wolstencroft & Smith |1978t Joos et al. 
2005). Thus, it seems that the tight anti-correlation between limb 
polarization and geometric albedo (within a small wavelength 
domain) is a generic feature of gas planets with a substantial 
amount of reflected light due to Rayleigh type scattering. 

According to simple Rayleigh scattering models a limb po- 
larization is expected for scattering atmospheres with some 
amount of Rayleigh-type scattering. 

However, for homogeneous (semi-infinite) Rayleigh scatter- 
ing atmospheres one does not expect an enhanced limb polar- 
ization in deep absorption bands (see Table 3 in paper 1) except 
for a small range of single scattering albedos around 0.8 and 1.0 
corresponding to geometric albedos between 0.3 to 0.79. A tight 
anti-correlation between polarization and albedo is obtained for 
models with finite Rayleigh scattering layers above diffusely re- 
flecting ground layers or deeper layers with different albedos 
(see Table 4 in paper 1). This can easily be understood as a con- 
stant amount of Rayleigh scattered radiation producing the limb 
polarization, combined with an albedo dependent contribution of 
unpolarized light from diffusely reflecting lower layers. 

The limb brightening effect points to a similar stratifica- 
tion in the atmosphere of Uranus and Neptune. As predicted 
by Belton & Price (U973I I and illustrated by e.g. Hammel et 
al. (1 19891 1 or Sromovsky (120051 1. only inhomogeneous scatter- 
ing atmospheres with a single scattering opacity increasing with 
depth will show a limb brightening. 

It is well known that Uranus and Neptune have optically 
thin scattering haze layers high in the atmosphere. This finding 
is based on careful analyses and modeling of the albedo spec- 
trum, center to limb profiles, and other multi-wavelength stud- 
ies of these planets (e.g. Baines & Bergstralh 1986, Hammel 
et al. [19891 Bain es & Smith [19901 Baines & Hammel [19941 
Sromovskv EOOBl 

With the newly detected limb polarization of Uranus and 
Neptune it would now be possible to test the existing mod- 
els with radiative transfer calculations including polarization. It 
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should not be expected that the atmospheric models have to be 
radically changed for Uranus and Neptune due to the polariza- 
tion signal. However, one may expect that an analysis including 
the scattering polarization would provide important new insights 
on the scattering properties for the aerosol particles in the haze 
layers of these planets. 
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